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SUMMARY 


An analysis is presented to calculate damping in ring seals for a compress- 
ible fluid. Results show that damping in tapered ring seals (optimized for 
stiffness) is less than that in straight bore ring seals for the same minimum 
clearance. Damping in ring seals can promote fractional frequency whirl and 
can, thus, be detrimental. Thus, tapered seals can benefit rotor and seal 
stability by having lower damping as well as higher stiffness. Use of incom- 
pressible results leads to large errors. 


INTRODUCTION 


Ring seals (annular seals) can have a considerable influence on the dynam- 
ic behavior of rotors. This is not surprising when one considers that such a 
seal has the appearance of a journal bearing, although with larger clearances 
than usual bearing practice. With a large pressure difference across the seal, 
most of the force generated is due to the high velocity throughflow of sealed 
fluid. In a series of papers (refs. 1 to 3) Black and coworkers calculated 
stiffness and damping in annular seals having constant clearance in the axial 
direction and sealing an incompressible fluid. In references 4 and 5, it was 
shown that stiffness could be considerably increased by tapering the seal bore 
so that clearance is greater at the inlet than at the exit. The higher stiff- 
ness can be beneficial in stabilizing rotors by shifting critical speeds; 
floating-ring seals can also benefit from higher stiffness with resulting 
longer life. Such a tapered bore configuration solved a wear problem in the 
hot gas seal of the Space Shuttle high pressure oxygen turbopump. 

Damping, as well as stiffness, influences seal and rotor behavior. As pre- 
viously mentioned, damping coefficients were calculated by Black, et al. 

(refs. 1 to 3) , for constant-clearance seals and incompressible fluids. Refer- 
ence 6 presented damping data for tapered ring seals with incompressible fluids. 
Heretofore no damping information has been available for seals passing compress- 
ible fluids. The purpose of this work is to provide such data for both 
constant-clearance and tapered seals. 
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SYMBOLS 


B damping coefficient 

- 3/3 

B dimensionless damping, BC^/RL \i Q 

B dimensionless damping, BL/Re^C^ 

C seal radial clearance 

c specific heat at constant volume 

v 

e seal eccentricity 

F seal force 

f Fanning friction factor 

g seal force per unit of circumference, eq. (40) 

H clearance ratio, 0^/C^ 

h seal local film thickness 

K seal stiffness 

K dimensionless stiffness, KC^/p^LD 

k entrance loss coefficient 

L seal length 

M Mach number, u/Vy^T 

p pressure ratio, p t /p g 

p pressure 

2 

Q M 

x s 

q heat flux 

R seal radius 

Re^ "sonic" Reynolds number, 2 Pq C 2\A^ T 0 /ta 0 
gas constant 

T absolute temperature 

t time 
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u 

velocity ratio, u /u 
t s 

u 

axial velocity 

V 

dh/<3t 

W 

Mach number ratio, M /M 

t s 

z 

axial coordinate 

a 

seal taper angle 

Y 

specific heat ratio 

e 

eccentricity ratio, e/C 2 


1 + k 

G 

circumferential coordinate 

M- 

dynamic viscosity 

P 

fluid density 

a 

density ratio, p /p 
t s 

9 

altitude angle, tan" 1 (Bcu^/K) 

O) 

P 

orbital frequency 

CD 

r 

shaft rotational speed 

Subscripts : 

e 

entrance 

rad 

radial 

s 

steady state 

t 

perturbed 

tan 

tangential 

X 

seal exit 

0 

upstream stagnation condition 

1 

seal entrance 

2 

seal exit 

3 

downstream reservoir condition 


ANALYS IS 


The configuration to be analyzed is that of a tapered bore ring seal whose 
clearance decreases in the flow direction (fig. 1). Around the circumference 
the clearance is given by 

h (z ,0) = C(z) + C 2 e cos 9 ( L ) 

where C is the clearance when the seal is concentric. The variation of C 
with axial coordinate z is 

C (z) = - az ( 2 ) 

The analysis is applicable to a straight seal by setting the taper angle 
a to 0. 

When the seal has a velocity relative to the shaft in the radial direc- 
tion, a damping force is generated. To determine this force, the pressure dis- 
tribution within the seal will be calculated and integrated over the seal area. 
The following assumptions are employed: 

(1) Eccentricity is small compared with the concentric clearance; that is, 
e «1. 


(2) The fluid behaves as a perfect gas. 

(3) Rotational effects on the flow field are neglected; the flow is one- 
dimensional in the axial direction. 

(4) The friction factor is constant everywhere within the seal. 

(5) Time derivatives of the fluid properties and velocity are neglected; 
that is the flow field is quasi-steady-state. Reference 6 examined this as- 
sumption for an incompressible fluid and concluded that acceptable accuracy re 
suits . 


The analysis begins from the time-dependent equations of continuity, 
momentum, and energy as presented by Shapiro (ref. 7). For a channel, which 
the seal passage approximates, these are, respectively, 

s <ouh> + & <ph) ■ 0 


q 



(pu 2 h) 



(4) 

(5) 


172 



The momentum equation may be simplified by using the continuity equation. In 
turn, the energy equation may be simplified through use of the continuity and 
momentum equations. Additionally, for a perfect gas 


and 


The seal taper angle is 


Defining 


£ = m 

p 



dz 


V 


__ dh 
~ St 


( 6 ) 

(7) 

( 8 ) 

(9) 


allows the continuity, momentum, and energy equations to be written as, respec- 
tively, 


a£ 

at 



+ U 


az 


+ 


p 


au 

az 



o 


( 10 ) 


q/h + 


ap = 

fpu 2 

du 

au 

az 

h 

+ p ^ + pu 

az 

fpu 3 £V _ _ 

1 

dg b£ 

VP 

h + h y 

- 1 

St dz 

P 


4- u 


( 11 ) 

( 12 ) 


The next step is to invoke the small perturbation assumptions and write 
the pressure, density, and velocity as the sum of a steady-state value plus a 
small perturbation 


P 

P 


P + P 


P„ + P, 


s t 7 

u = u 4- u 
s t 


The steady-state equations are 

d P. 


du 


s , s 

u ~a — + p ~r~ 

s dz s dz 


p u a 

s s 


dz 


fp u 

s s 

h 


4- p u 
s s 


du 
_ s 

dz 


(13) 


(14) 


(15) 
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, 3 

fp u 

/, s s 

q s /h + “h - 


- 1 \ dz 


d P s YP S dp s 


The steady-state equations were solved (for q = 0) in reference 5 and will not 
be dealt with further. The perturbed equations of continuity, momentum, and 
energy are, respectively, 


t s , d , , x 

: — + x- (p u + p u ) - 
St h dz VH t s s t 


p u a p u a 

St t s 


Sp fp u 2 fp u u 
_ H t s s s t 

Sz h h 


ou t au t ou 

+ P 37* + p u 3 — + (p u + P u ) 

s dt s s dz s t t s o* 


3 2 

fp . u 3fp u u p v 
„ K t s s t s , s 

q /h + — - — + r + —r— 

M t h h h 


1 K 5p s YP t ^ s *P 8 ^ s P t U s 

- 1 dt + U s dz + U t dz p dz + 2 dz 

L s p s 


YP S /dP t 3 

p \ dt + U s dx + U t dz 


In accordance with the assumptions, the time derivatives of p t , p t , and 
u t are neglected, and the heat flow q is taken to be zero. Equations (17) 
to (19) then become ordinary differential equations which may conveniently be 
written in dimensionless variables 


p-^ 

p s 


Additionally, the Mach number is defined as 


As with the other variables, the Mach number may be considered as the sum of a 
steady-state quantity and a small perturbation. 
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(22) 


Defining 


we find, from equation (21), 


M = M + M 
s t 




(23) 


(24) 


After some algebraic manipulation, making use of equations (14) to (16) , the 
perturbed continuity, momentum, and energy equations become, respectively, 


wherein 


da 

dz 


dU v 

dz u h 
s 


0 


dP dU = 2W ^s 

dz ^ dz p dz 

*• o 


dP 

dz 



■ ^ (w + Jr) 


Q 



(25) 


(26) 

(27) 


(28) 


Thus, we have three simultaneous, ordinary differential equations in four 
unknown variables P, a, U, and W. A fourth equation needed to obtain a solu- 
tion may be obtained by differentiating equation (24) with respect to z: 


dW _ dU 1 /do dP\ 

dz dz + 2 \dz ~ dz/ ( 29 ) 


Equations (25) to (27) and (29) may now be solved simultaneously for the 
perturbed variables. As for the steady-state variables (ref. 7), the differen- 
tial equations for the perturbed variables may be combined to yield a single 
differential equation in the perturbed Mach number W. 

Reference 8 derives a differential equation for the steady-state pressure 
p s , which may be written as 



h('l Y - Q) Ca + f + f (Y - 1)Q] 


(30) 


With this, the equation for perturbed Mach number is 
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dW 

dz 


W 

h(l - Q) 


( (Y Z ^ [a + f + f(Y - 1)Q] 


+ (y - DfQ(vQ + l) } 


(Y + l)v 

2yuh(l - Q) 


(31) 


From the 
Equation (31) 
eq. (2)) 


steady-state solution (ref. 5), Q is a known function of 
will be easier to solve if we make the substitution (ref. 


z . 

5, 


^ = 
dz 



Y 2 ~ q)(y £ Q + ») 
h(l - Q) 


(32) 


Also, u s varies according to (ref. 8) 



(33) 


where the subscript e denotes the steady-state value at z = 0. Equation (31) 
becomes 


dW 

dQ 


W 

2(1 + Y ~ 1 q)(y £ Q + a) 


Y .+ J: 
1 - Q 


a + 



l ± vQ 

1 - Q 


+ YQ + 1 



(Y + l)vM e 

4yu e Q 3/2 (yfQ + a) |/ (l + q) (l + 


(34) 


The perturbed Mach number W may now be found by solving equation (34). 

To determine the seal damping force, the perturbed pressure P must be 
known. Equations (25) to (27) may be combined to yield a differential equa- 
tion for P. Making use also of equations (32) and (33) yields 


dP 

dQ 


yW 


1 + 


l) (YfQ 


+ a) 


a I + . Y Q - Q + f(l + YQ - 



vM 

+ 6 — (35) 

2u e M s (YfQ + a) fA 1 + Q ) l 1 + Q e) 

Boundary conditions . - At the seal exit it is assumed that there is no 
change in the boundary conditions of the steady-state problem. Thus, for 
choked flow 
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For flow which is not choked 



(36) 


(37) 


At the seal entrance, the pressure and Mach number are related by (ref. 5) 


JL = 


P 


0 



\Y/(1-y) 

IQ 


(38) 


The pressure and Mach number may be written as the sums of their steady-state 
and perturbed components and a binomial expansion performed on the right side 
of equation (38). After neglecting all powers of Q higher than unity and 
subtracting out the steady-state terms, there remains 


P = - 


YhQgW 


1 + 


(39) 


. Solution of equations. - The results presented herein were obtained 
through a numerical solution of the differential equations for W and P 
(eqs . (34) and (35)). For the case of choked flow, Q = 1 at the seal exit- it 
is not possible to integrate numerically to this limit because of the term’ 

1 - Q m the denominator of equations (34) and (35). However, W and P may 
be determined as accurately as desired by taking the solution close to Q = 1 
without reaching the limit. To ascertain the validity of this approach, equa- 
tions (34) and (35) may be rewritten, in simplified form, for the case when 
Q - 1. The resulting equations may be solved analytically, yielding solutions 
° r ^ w hich remain finite as Q -+ 1. Thus 3 the error is of the order 

( ^lim'^Q = 1 which can be made as small as desired by taking Oi • arbi- 
trarily close to 1. ViLm 

. Cg L . cu l- a tj; Qn of s e al forces . - The seal force per unit of circumference due 
to the perturbed pressure p t is given by 


g 




(40) 


where Q x is the value of Q at z = 
not contribute to the total seal force 
ference. From reference 8 


L. The steady-state pressure p s does 
because it is uniform around the circum- 
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Making use of this and equation (32) results in 

Q, 


S = lP e h e M er L+ 2 


V 1 + ^ / 


'X 


pq - Q) d Q 


Q 3/2( l+ii J p ) 3/2 (v£Q 


a) 


^^0 integral in equation (42) is easily evaluated numerically. 
The total seal damping force in the direction 9=0 is 


(41) 


(42) 


/ 2 it 

gR cos 0 d 0 

0 


For a seal velocity in the direction 9-0 


and , thus , 


v = Vq cos 9 


g = g 0 cos e 


where the subscript 0 now denotes the condition at 0 
integral in equation (43) may be evaluated, yielding 


F = -TrRg 


0 


A damping coefficient is now defined: 


F ^0 


(43) 


(44) 

(45) 


0. With this the 


(46) 


(47) 


or 


B 


phM 
2 Vq e e e 




pq - Q) d Q 

3/2 

1 2 ~ q) (vfQ + a ) 


(48) 


where it is understood that P is evaluated for v = v Q . The actual value of 
v 0 is immaterial. Since both W and P are directly proportional to v 0 , 
it will cancel out of the expressions for the proper choice of dimensionless 
variables . 
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RESULTS 


The analysis of the 
puter. The differential 
integrator with automatic 
results of reference 5. 


preceding section was implemented on a digital corn- 
equations were solved using a fourth-order Runge-Kutta 
error control. The computer programs made use of the 


Results are presented in figures 2 and 3 for choked flow and in figures 4 
an or a sealed pressure ratio of 2, for which the flow is not choked. The 
independent variable in all c ases is Re 0 C 2 /L where Re 0 is the -sonic" 
Reynolds number 2C 2 p 0 V7 *t£/h 0 . This independent variable was chosen be- 
cause laminar flow results are then independent of seal clearance- to- length 
ratio C 2 /L. The specific heat ratio y was 1.4 for all cases. Results are 

shown for straight seals and for tapered seals optimized for maximum stiffness- 
to-leakage ratio (ref. 5). 

d o /T gUr ^ 2 shows laminar How results for small values of seal parameter 
ReQC 2 /L. Damping is plotted in terms of the dimensionless variable 


bc; 




(49) 


There is little variation for the range of seal parameter up to 10. Damping 
approaches a constant value as Re 0 C 2 /C - 0. Damping is higher in the 
straight seal than in the tapered design because of the lower average clear- 
ance (seals are compared on the basis of minimum clearance). 

Figure 3 shows results for choked flow over the complete range of seal 
parameter for both laminar and turbulent flow. In this figure, damping is 
shown in terms of the dimensionless quantity 


B = 


BL 


BC, 


Re 0 C 2 2RL 2 V^YPqPq 


(50) 


For this choice of variable, laminar flow damping varies strongly with seal 
parameter. In contrast, turbulent flow damping shows a much smaller variation 
For turbulent flow, damping depends on clearance-to-length ratio Co/L, appear- 
ing to reach a minimum for the middle values of C 2 /L investigated, and rising 

for either higher or lower values. Overall, however, variation is small even 

for the extreme values of C 2 /L considered; thus, for all other factors held 

constant, equation (50) indicates that damping is approximately proportional to 
the square of seal length. 

The left end of the curves for turbulent flow corresponds to a Reynolds 
number m the seal passage of 3000; this is generally considered the lowest 
^ e . for which one can be assured of turbulent flow. A Reynolds number of 
2300 is usually taken as the upper limit for laminar flow. Points where 
Re - 2300 are shown for various C 2 /L values on the laminar flow curves. 
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Damping for a pressure ratio of 2 is shown in figure 4 for small values of 
Re„C,/L For this pressure ratio the flow is not choked. In common with fig- 
ure 2, damping does not vary greatly with seal parameter. 

Results for a pressure ratio of 2 are shown in figure 5 for the complete 
range of seal parameter. As in figure 3, the damping is expressed as 


BC, 


B = 


2 

2RL VSt 0 P 0 


For the straight seal, damping values are lower than for choked flow, as less 
mass passes through the seal when the flow is unchoked. For the tapered seal, 
damping for p 0 /p 3 = 2 is higher than the choked case when the flow is lam ^ ar 
and somewhat l2we? when the flow is turbulent. This apparently anomalous be- 
havior for laminar flow can be explained by observing that the optimum clear- 
ance ratio is lower for p 0 /p 3 = 2 than for choked flow; thus, the average 
clearance is less than in the seal optimized for choked flow. Optimum clear- 
ance ratios H = C]/C 2 appear in figures 6 and 7 (taken from ref. 5). 

Example of seal calculation . - Calculate damping for the example seal of 
reference 5. Pertinent seal data are 


L = 10 mm 

D = 50 mm 

C 2 = 0.05 mm 

p Q = 100 bar 

P 3 = 4 bar 

T q = 800 K 

Fluid: Hydrogen gas 


Thus, it was calculated in reference 5 that 

^2 

— = 0.005 

L 


The flow is choked; 
seal 

yielding 


thus , 


Re 0 = 

Re ^ 

R 0 L 


damping data 


37 100 
= 186 

are taken from figure 


3. 


B = 0u89 


B = 570 N sec/m 


For a straight 
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For the tapered seal 


B = 0.47 

and 


B = 300 N sec/m 


Effect of seal rotation . - The effect of rotation was examined in refer- 
ence 6 for a seal passing incompressible fluid using an approximation derived 
y Black and Jenssen (ref. 2). For the analysis presented herein, fluid tem- 
poral derivatives were neglected; thus, the inertia terms calculated in refer- 

ence 6 do not appear in the present results. The expression for seal forces 
then becomes 



F 


y 


K 




In equation (51) CD,, is the shaft rotational speed. 


(51) 


In the usual case of seal motion, the seal journal describes some orbit 
wi thm the seal clearance. We will consider the case of a centered circular 
orbit with orbital frequency o>p. For this condition, it is easy to calculate 
the radial restoring force and the tangential (whirl direction) force. Similar 
to reference 6, we find 


rad 


= -Ke 


tan 


= Be 


a 


to - to 


(52) 


where e is the eccentricity of the seal (fig. 8). The radial force depends 
only on the eccentricity e (the minus sign appears because F rad is defined 
as positive in the outward direction). The tangential force, however, depends 
on the relative magnitudes of the spin and whirl speeds. Equation (52) may be 
written in terms of radial and tangential stiffnesses: 


K 


tan 


K 


rad 






(53) 


A positive value of K tan will inhibit positive whirl; conversely, negative 
K tan promotes positive whirl. Equation (53) for K tan shows that a positive 
damping coefficient promotes forward whirl when aip/a* < 1/2 and inhibits whirl 
when cnp/cu,- > 1/2. As pointed out in reference 6, this is remarkably similar 
to the behavior of a full circular journal bearing. The seal differs from a 
self-acting bearing (but is similar to an externally pressurized bearing) in 
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that the radial stiffness is independent of whirl speed. Equation (53) also 
shows that damping in a rotating seal is not entirely beneficial. Damping pro- 
duces a tangential force which inhibits synchronous whirl but promotes whirl at 
half frequency or less. 

Example of whirl forces. - Seal damping has been presented in dimension- 
less form, as were stiffness and leakage in reference 6. Also, damping force 
varies with rotor speed (eq. (53)). Thus, the relative size of the radial and 
tangential stiffness of equation (53) can only be compared for specific exam- 
ples. The seal example above is similar to one of the hot-gas seals in the 
Space Shuttle high-pressure oxygen turbopump. This pump has a nominal operat- 
ing speed of 31 000 rpm at full-power level. For this speed, table I compares 
the stiffness and damping for straight and tapered seals. Stiffness data are 
taken from reference 6. As table I shows, tapered seal stiffness is much high- 
er than that for a straight seal, and damping lower. This means that for a 
tapered seal the force generated is more nearly in line with the displacement, 
a situation that is considered to promote stability. 

Difference between compressible and incompressible results . - Before the 
availability of stiffness, damping, and leakage flow data for seals handling 
compressible fluids, some workers used incompressible data as the "best avail- 
able." It is instructive to examine the difference in results using compress- 
ible and incompressible analyses. Table II compares results for the seal used 
in the example above. Compressible fluid stiffness and leakage flow are taken 
from reference 5. Incompressible data from references 4 and 6 are applied two 
ways: first using the upstream (stagnation) fluid density and second using 

fluid density as the mean of upstream and downstream densities. 

The incompressible theory grossly overestimates seal stiffness, by an 
order of magnitude for the straight seal. There is little difference in using 
upstream or mean density. Damping is underestimated by the incompressible the 
ory more severely when the mean density is used. Leakage flow is overesti- 
mated by incompressible theory, although not as much when a mean density is 

used . 


One would not expect incompressible theory to yield accurate results be- 
cause the nature of the flow is much different than with a compressible fluid. 
Fluid compressibility results in decreasing density and increasing velocity in 
the flow direction, and the flow becomes choked for large pressure ratios. The 
pressure gradient increases in the flow direction and, theoretically, becomes 
infinite at the seal exit for choked flow. 


CONCLUDING REMARKS 


An analysis has been performed to calculate the damping in straight and 
tapered seals for a compressible fluid. Results show that damping in optimized 
tapered seals is considerably less than in straight seals for the same minimum 
clearance. It was also pointed out that damping in rotating seals is not en- 
tirely beneficial as it can sometimes promote whirl. Thus, for rotor or seal 
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stabilization, tapered seals may confer a double benefit. Not only is tapered 
seal stiffness generally higher than that for straight seals, but damping is 
lower. 
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TABLE I. - STIFFNESS AND DAMPING FOR HPOTP 


STRAIGHT AND TAPERED SEALS 



Straight 

seal 

Tapered 

seal 

Dimensionless stiffness, K 

0.015 

0.095 

Stiffness, K, MN/m 

1.5 

9.5 

Dimensionless damping, B ! 

0.89 

0.47 

Damping, B, N sec/m 

570 

300 

BCD , MN/m 
r 

1.9 

1.0 

Attitude angle <p 5 deg 

51 

6 


TABLE II. - COMPARISON OF COMPRESSIBLE AND INCOMPRESSIBLE RESULTS FOR 


HPOTP STRAIGHT AND TAPERED SEAIS 



Stiffness 

, K, MN/m 

Damping, B 

, N sec/m 

“Leakage flow, g/sec 


Straight 

seal 

Tapered 

seal 

S traight 
seal 

Tapered 

seal 

S traight 
seal 

Tapered 

seal 

Compressible theory 

1.5 

9.5 

570 

300 

17 

22 

Incompressible theory, 
stagnation density 

15.1 

29.1 

308 

245 

32 

42 

Incompressible theory, 
mean density 

14.6 

28.6 

223 

176 

22 

29 







OPTIMUM TAPER 

STRAIGHT 


O Re « 2300 FOR GIVEN C?/L 



SEAL PARAMETER, Re 0 (C 2 /LI 


Figure 3. - Damping in ring seals; choked flow, laminar and turbulent. 



SEAL PARAMETER, Re 0 (C 2 /L) 

Figure 4. - Damping in ring seals; unchoked, P 0 /P 3 = 2, 
laminar flow, small values of seal parameter. 


OPTIMUM TAPER 

STRAIGHT 

\ O Re ■ 2300 FOR GIVEN C 2 /L 



SEAL PARAMETER. Re 0 <C 2 /L) 

Figure 5. - Damping in ring seals; pg/p^ = 2 , laminar and turbulent flow. 
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Figure 6. - Optimum clearance ratio of tapered seals, choked flow. 
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Figure 7. * Optimum clearance ratio of tap< 
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